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Quantum vacuum forces dictate the interaction between individual atoms and dielectric 
surfaces at nanoscale distances. For example, their large strengths typically overwhelm 
externally applied forces, which makes it challenging to controllably interface cold atoms with 
nearby nanophotonic systems. Here we theoretically show that it is possible to tailor the 
vacuum forces themselves to provide strong trapping potentials. Our proposed trapping 
scheme takes advantage of the attractive ground-state potential and adiabatic dressing with 
an excited state whose potential is engineered to be resonantly enhanced and repulsive. This 
procedure yields a strong metastable trap, with the fraction of excited-state population 
scaling inversely with the quality factor of the resonance of the dielectric structure. We 
analyse realistic limitations to the trap lifetime and discuss possible applications that might 
emerge from the large trap depths and nanoscale confinement. 



1 ICFO— Institut de Ciencies Fotoniques, Mediterranean Technology Park, Castelldefels, 08860 Barcelona, Spain. 2 Joint Quantum Institute, College Park, 
Maryland 20742, USA. 3 National Institute of Standards and Technology, 100 Bureau Dr MS 8410, Gaithersburg, Maryland 20899, USA. 4 IQIM, California 
Institute of Technology, Pasadena, California 91125, USA. 5 Norman Bridge Laboratory of Physics 12-33, California Institute of Technology, Pasadena, California 
91125, USA. Correspondence and requests for materials should be addressed to D.E.C. (email: darrick.chang@icfo.es). 



NATURE COMMUNICATIONS | 5:4343 | DOI: 10.1038/ncomms5343 | www.nature.com/naturecommunications 

© 2014 Macmillan Publishers Limited. All rights reserved. 



1 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms5343 



One of the spectacular predictions of quantum electro- 
dynamics is the emergence of forces that arise purely 
from quantum fluctuations of the electromagnetic 
vacuum 1 . Also known as London-van der Waals 2 or Casimir 
forces 3 in different regimes, these forces are often dominant at 
short distances and can give rise to undesirable effects in 
nanoscale systems, such as nanomechanical stiction 4 . These 
forces have also attracted increasing attention in the fields of 
atomic physics and quantum optics. In particular, significant 
efforts have been made in recent years to interface cold atoms 
with the evanescent fields of dielectric micro- and nano-photonic 
systems 5-14 . These systems are expected to facilitate strong, 
tunable interactions between individual atoms and photons for 
applications such as quantum information processing 15 and the 
investigation of quantum many-body physics 16-18 . In practice, 
efficient atomic coupling to the evanescent fields of the 
nanophotonic systems requires that atoms be trapped within 
sub-wavelength distances of these structures. 

At these scales, quantum vacuum forces can overwhelm the 
forces associated with conventional optical dipole traps, typically 
resulting in a loss of trap stability at distances d< 100 nm from 
dielectric surfaces. Given the ability to engineer the properties of 
nanophotonic structures, an interesting question arises as to 
whether such systems could be used to significantly modify 
vacuum forces, perhaps changing their sign from being attractive 
to repulsive, or even creating local potential minima. The strength 
of nanoscale vacuum forces should lead to unprecedented energy 
and length scales for atomic traps, which would find use beyond 
nanophotonic interfaces, such as in quantum simulation proto- 
cols based on ultracold atoms 19 and control of inter- atomic 
interactions . 

Here we propose a novel mechanism in which engineered 
vacuum forces can enable the formation of a nanoscale atomic 
trap. While a recent no-go theorem 21 forbids a vacuum trap for 
atoms in their electronic ground states, tailored nanophotonic 
systems can yield strong repulsive potentials for atomic excited 
states. We show that a weak external optical field can give rise to 
an overall trapping potential for a dressed state. Remarkably, 
absent fundamental limits on the losses of the surrounding 
dielectric structure, the fraction of excited- state population in the 
dressed state can become infinitesimal, which greatly enhances 
the trapping lifetime and stability. We identify and carefully 
analyse the actual limiting mechanisms. While we present 
calculations on a simple model where analytical results can be 
obtained, we also discuss the generality of our protocol to realistic 
systems such as photonic crystal structures. 

Results 

Quantum interactions between atom and surface. The no-go 
theorem for non-magnetic media states that a dielectric object in 
vacuum cannot be stably trapped with vacuum forces for any 
surrounding configuration of dielectric objects, provided that the 
system is in thermal equilibrium 21 . In analogy with Earnshaw's 
theorem, which prohibits trapping of charged objects with static 
electric potentials, at best one can create a saddle-point potential 
(see, for example, ref. 22 involving metal particles and refs 13,23 
for a hybrid vacuum/optical trap for atoms). For dielectric 
objects, potential loopholes involve embedding the system in a 
high-index fluid 24 or using blackbody radiation pressure 
associated with strong temperature gradients 25 . Applied to 
atoms, this theorem essentially forbids the stable vacuum 
trapping near a dielectric structure when the atom is in its 
electronic ground state. No such constraint exists for atoms in 
their excited states 26 , although the robustness of any possible trap 
would be limited by the excited- state lifetime. Here, we show that 



an atom only weakly dressed by its excited state can be trapped by 
properly tailoring the dispersion of the underlying structure. 

We begin by briefly reviewing how quantum vacuum 
fluctuations give rise to forces on the ground and excited states. 
Within an effective two-level approximation of an isotropic 
atom with ground and excited states |g>, |e>, respectively, the 
interaction between the electromagnetic field and an atom at 
position r is given by the dipole Hamiltonian H— — d-E(r). 
Conceptually, following a complete decomposition of the 
electric field into its normal modes k, one can write 

H = - 2k&( r )(£eg + £ g e)(«k + flk)> where °ij = \ i ) ( J\- The 
energy non- conserving terms (<J eg al an d 0" ge 0k) enable an atom 
in its ground state \g, 0 ) to couple virtually to the excited state 
and create a photon, \e, lk), which is subsequently re-absorbed. 
The corresponding frequency shift for the ground state within 
second-order perturbation theory is <5co g (r) = — J] k gk(r) 2 / 
(cL>o + ^k)> where co 0 is the unperturbed atomic transition 
frequency. This shift results in a vacuum- induced mechanical 
potential when translational symmetry is broken owing to 
dielectric surfaces. Using a quantization technique for electro- 
dynamics in the presence of dispersive dielectric media, the shift 
can be expressed in terms of the scattered component of the 
dyadic electromagnetic Green's function evaluated at imaginary 
frequencies co = iu (ref. 27) (also see Methods), 

oo 2 

^ W =3 of / dM cofT^ TrGsc(r ' r ' /M) ' (1) 

0 

where T 0 is the free-space spontaneous emission rate of the atom. 
For simplicity, we have given the result in the zero-temperature 
limit. Finite-temperature corrections 28 ' 29 associated with atomic 
interactions with thermal photons are only relevant for distances 
from surfaces comparable with or larger than the thermal 
blackbody wavelength, d>A T , where = yy ~ 7.6 \im at 
T=300K. 

A similar process can occur for an atom in its excited state 
\e, 0), in which the atom virtually emits and re-absorbs an off- 
resonant photon (thus coupling to the state \g, lk))- However, 
the excited state is unique in that it can also emit a resonant 
photon. The total excited-state shift is given by 27 

S(D e (r) = - -^TrReG sc (r,r,co 0 ). (2) 

3 ojq 

The first term arises from off-resonant processes, while the 
second term can be interpreted as the interaction between the atom 
and its own resonantly emitted photon. Naturally, the dielectric 
environment can modify the spontaneous emission rate as well, 

r(r) = T 0 + 2r ° nc Tr Im G sc (r, r, co 0 ) . (3) 

Note that the resonant shifts and modified emission rates are 
complementary, in that they emerge from the real and imaginary 
parts of the Green's function. Micro- and nano-photonic systems 
are already widely exploited to modify emission rates of atoms 
and molecules in a variety of contexts 30-33 . We propose that these 
same systems can be used to tailor excited- state potentials to 
facilitate vacuum trapping. 

Interactions with a Drude material. In what follows, we present 
a simple model system that enables one-dimensional trapping in 
a plane parallel to a semi-infinite dielectric slab, as illustrated in 
Fig. 1. The use of a homogeneous dielectric (as opposed to, for 
example, a photonic crystal) greatly simplifies the calculation 
and enables one to understand the relevant trap properties 
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Figure 1 | Schematic of vacuum trapping mechanism, (a) A semi-infinite dielectric with specially tailored frequency-dependent permittivity e(co) enables 
large repulsion in the excited state co e (z), in contrast to the attractive vacuum potential experienced by an atom in its ground state. An overall 
trapping potential normal to the surface (along z) can be generated by driving the atom with an external laser of frequency a> L , which is greater than the 
natural resonance frequency co 0 of the atom. The laser comes into local resonance with the atom at z b . An atom far from z b is unaffected owing to the large 
laser detuning co L -co 0 and feels an attractive force F toward the dielectric interface, (b) The atom comes closer to local resonance with the laser 
as z->z b . This creates an atomic-dressed state, whose excited-state component yields a net repulsive force on the atom, (c) A slowly moving atom 
experiences the adiabatic dressed-state potential U d (z) shown in red, which characterizes the position-dependent mixing between ground and 
excited states. The minimum of the dressed potential is located at z t >z b . 



analytically, although we later argue why the qualitative features 
should still hold in more complex, realistic settings. In the 
dielectric slab model, dispersion engineering will be realized via 
the frequency-dependent electric permittivity e(co). While this 
particular system achieves a one-dimensional trap along z, our 
electrodynamic calculations are performed in three dimensions. 

A short-distance expansion of the Green's function, valid for 
sub-wavelength scales, shows that the excited-state emission rate 
T(z) and resonant contribution Sco^ (z) to the shift are given by 



c5a4 r) (z) 



1 



4(M 



rim 



r„ 



8(M 



Re 



£a+l 



(4) 



(5) 



where the permittivity e a — e(co 0 ) is evaluated at the atomic 
resonance frequency, z is the distance from the surface and 
k 0 = co 0 /c is the resonant free-space wavevector. A large repulsive 
potential is generated when e a approaches — 1 from above, 
e a -> — 1 + , which is analogous to the interaction between a 
classical dipole at position z and its large induced image dipole in 
the dielectric. The modified spontaneous emission rate originates 
from the absorption or quenching of the atomic emission owing 
to material losses. 



We choose a Drude model, e(co) = 1 



as a simple 



co 2 + reoy' 

dielectric function. Note that this function satisfies Kramer s- 
Kronig relations (causality), as a choice of a non- causal function 
could result in an apparent violation of the no-go theorem 21 as 
well. In the limit of vanishing material loss parameter y, the 
system passes through e = — 1 at the plasmon resonance 
frequency co p \ = co p /\^2. The system is conveniently 
parameterized by the quality factor Q=(D p \ly and a 
dimensionless detuning A p = (co — (o p i)/y. For Q»A p , the 
factor ~ — Q/(/ + 2A p ) resembles the complex 

susceptibility of a simple resonator, which in the far- detuned 



i Q an d mi 

16(fc 0 z) A p r ° 



16(fc 0 z) 3 A p 



£. The 



limit (A p »l) yields ^ 

dispersion and dissipation scale like A" 1 and A~ 2 , respectively. 
Thus, for high Q, it is possible to choose detunings where the 
atom still sees significant repulsive forces, but where spontaneous 
emission into the material is not significantly enhanced (in 
addition to material-induced emission of equation (4), there still 
is emission into free-space at a rate ~r 0 ). This scaling behaviour 
is in close analogy with conventional optical trapping of atoms, 
where simultaneously using large trapping intensity and detuning 
maintains reasonable trap depths but suppresses unwanted 



photon scattering. The ground- state shift for atomic frequencies 
(D 0 ~(o p i and in the near-field is given by Sco g ^ — 32 ^°^3 . An 
additional condition for the validity of the short- distance results 
obtained for the excited state is k 0 z < ^/A p /Q, as shown in the 
Methods. We continue to present approximate results in this 
regime owing to its simplicity, although all of our numerical 
results (such as in Figs 2c,d and 3) are calculated using the full 
Green's functions. 

We also briefly comment on two effects in realistic systems that 
can be safely excluded from our analysis. First, in principle, the 
excited- state shift can be affected by coupling to even higher 
electronic levels. However, a careful calculation of this contribu- 
tion 34 shows that it is negligible compared with the resonantly 
enhanced term of equation (5), provided that we work in the 
regime Q » A p . Second, we have ignored corrections to the shifts 
owing to surface roughness. It can be shown that the length scale 
for modifications to the shifts is roughly set by the size of the 
surface variation itself 35 . The realistic regimes of operation of our 
trap will be at d> 10 nm from surfaces, which allows us to safely 
ignore the angstrom-scale surface roughness achievable in state- 
of-the-art dielectric structures 36 . 

Formation of a vacuum trap. We now describe how a trap can 

form for an atom subject to these potentials, and in the presence 

of a weak driving laser of frequency co L and Rabi frequency Q (see 

Fig. 1). For conceptual simplicity, we choose the Rabi frequency 

to be spatially uniform, such that all atomic forces are attributable 

to the vacuum potentials alone. Qualitatively, the strong position 

dependence of the ground- and excited- state frequencies causes 

the laser to come into resonance with the atom at a single, tunable 

point z — z^. Under certain conditions (described in detail later), 

an atom starting at position z » Zb will be far detuned, such that 

it is essentially in the ground state and is attracted by the pure 

ground- state potential. However, moving closer to Z\> brings the 

atom closer to resonance. The dressing of the atom with a small 

fraction of excited- state population causes a repulsive barrier to 

form near z~z b , yielding a metastable trap. Significantly, our 

approach does not attempt to directly counteract the attractive 

ground- state potential with large external optical potentials, in 

sharp contrast with other trapping schemes near dielectric sur- 
faces 8,9,ll,12,37-39_ 

The order of our calculation of the trap properties is as follows. 
We first find the adiabatic potential experienced by a slowly 
moving atom. We derive relevant properties around the location 
of the trap minimum z t (Fig. 1). In particular, we show that 
absent material constraints (arbitrarily high Q), the dressed state 
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Figure 2 | Trap properties, (a) Ideal limit of trapping potential, corresponding to the pure ground-state potential ftdco g (z) and an infinite barrier at z = z b . 
The maximum depth of the potential is denoted by U max =fi|^ft)g(z b )|. The quantum wavefunction amplitudes |i/Kz)| 2 (in arbitrary units) for the two 
lowest eigenstates are illustrated in red, while the binding energy U 0 and position uncertainty Az of the ground state are depicted as well, (b) Properties of 
the ideal trapping potential versus barrier position z b , for the case of a Caesium atom. The blue solid and dashed curves give the potential depth 
at the barrier position, U max , and the ground-state-binding energy U 0l respectively. The green curve depicts the ground-state uncertainty Az. (c) For the 
non-ideal trapping potential, the depth of the classical potential Ud ep th is numerically calculated. Here we plot it normalized by the maximum possible 
value U max =S|^ft)g(z b )|, versus detuning parameter A p and normalized Rabi freuqency Q = Q/Q m]n . (d) The trapping potential, ground-state-binding 
energy U 0l and wavefunction amplitude \\j/(z)\ 2 for the parameter sets A, B, as illustrated in c. The colour coding is the same as in a. 
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Figure 3 | Trap lifetime, (a) Trap lifetime T to tai versus dimensionless Rabi frequency Q = Q/Q m j n and detuning A p . Here we have assumed a 
plasmon resonance quality factor of Q = 10 7 and atomic properties corresponding to Caesium. The regions of parameter space where no trap exists are 
labelled 'NT' in the plot, (b) For the same parameters, the individual contributions of recoil heating, anti-damping and transient heating are plotted 
as a fraction of the total heating rate. Tunnelling contributes less than 10 ~ 2 to the total heating over this range of parameters and is not plotted here. 



can in principle have infinitesimal excited- state population and 
scattering. We then quantize the motion to find the motional 
eigenstates, binding energy and position uncertainty Az. Finally, 
we analyse the mechanisms that limit the trap lifetime. 

Following standard procedures 40 , the average force 
experienced by the atom is given by the expectation value of 

the Heisenberg equation ^ = — h (^r&ee(t) + ^^gg(0)> while 
the internal atomic dynamics satisfies the usual Bloch equations, 
for example, 

^p= (^a(^)-r(z)/2)0- ge + /y(0- e e-0- gg ). (6) 

Here, for notational simplicity, = (o^) denotes the expecta- 
tion value of the atomic variables, and we have defined S^z) = 
a> L — (a> e (z) — (i) g (z)) as the detuning between the laser and the 
local atomic resonance frequency. When the atom moves 
slowly on the timescales of the internal dynamics, the atomic 
coherence can be adiabatically eliminated, ^ « 0, such that 
the atomic populations are functions of position alone, 
cr ee — ft^2 — A 9, ,2 — We are primarily interested in the regime 

T(z) +4c> a (z) + 212 r ' ° 

where the atom is weakly driven, cr ee «l. A more careful 



calculation reveals that the Green's function G sc (r, r, co) for the 
excited-state resonant shift and emission rate in equations (2) and 
(3) must in fact be evaluated at the laser frequency in this regime, 
which reflects that the atom primarily acts as a Rayleigh scatterer 
(see Methods). Consequently, e a in equations (4) and (5) is 
replaced with e(co L ). The adiabatic potential seen by a slowly 
moving atom is U&(z) — dz d ^\ In principle, the Heisen- 
berg equation for ^ should also include the effects of spontaneous 
emission. It can be shown, however, that this term contributes a 
zero average force (see Methods), and only results in a heating 
term to be analysed later. 

Three independent frequencies characterize the system (co L , co p 
and cd 0 ), and two relative frequencies must be specified. One is 
determined by selecting the trap barrier position, determined by 
the relation (5 a (z b )=0. A second is the dimensionless detuning 
between the plasmon resonance and laser, A p = (a> L — co p i)/y, 
which we treat as a free parameter that determines the relative 
strengths of the excited-state shift and dissipation. In the near- 
field, the excited- state repulsion exceeds the ground- state 
attraction by a factor 2Q/3A p . The population at z h must then 
exceed the inverse of this quantity, <7 ee >3A p /2Q, to provide a 
classical barrier in the adiabatic potential, which translates to a 
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minimum Rabi frequency « 1 , /7 1 , 3 * Since A n < Q, this 

^ ; r 0 i6(fc 0 Zb) V p 

scheme potentially enables atom trapping near surfaces with 
greatly reduced intensities as compared with conventional 
trapping. In the latter case, an optical dipole potential requires 
a minimum driving field of ~ „ 1 , 3 to overcome the 

° 1 o (koz h ) 

attractive ground- state potential. 

The excited-state population at the trap minimum is fixed by 
the ratio of ground to excited-state forces, cr ee (z t ) = and is 
surprisingly independent of Q (provided that Q>Q min ). An 
increasing Rabi frequency has the effect of increasing the local 
detuning S a (z t ) to preserve the same population. This manifests 
itself as a self- selection of the trapping position and an increasing 
distance from the barrier with increasing intensity, 



Zt-z\> 



1 

6Av 



Q 2 -l 



1/2 



(7) 



where we have defined the dimensionless Rabi frequency 
Q = Q/Q min . The ability to tune the trap position using both 
the laser amplitude and its frequency (which determines z b ) 
enables the adiabatic transformation between deep but small traps 
close to the surface to larger, shallow traps ~100nm from the 
surface. This is expected to greatly facilitate trap loading, as 
techniques to trap atoms at these larger distances are already well 
established 8 ' 11 ' 12 ' 41 . 

The photon-scattering rate at the trap position, R sc = 
r(z t )cr ee (z t ), leads to lifetime limitations of our trapping scheme. 
Fixing other parameters, the detuning A p from the plasmon 
resonance can be optimized to yield the minimum scattering rate, 



R. 



Thus, absent fundamental limits on the 



~ 4V2Q(te) 3/2 ' 

magnitude of Q, the photon -scattering rate can in principle be 
suppressed to an arbitrary degree. This implies that an 
infinitesimal violation of the assumptions underlying the no-go 
theorem can in fact allow for trapping. 

In a conventional optical dipole trap, the potential U— 
— a(co L )|£(r, cl> l )| 2 /2 seen by an atom depends on the local field 
intensity and is proportional to a spatially constant atomic 
polarizability evaluated at the laser frequency. In contrast, here we 
exploit a novel mechanism arising from strong spatial variations 
in the polarizability itself, induced through the vacuum level 
shifts. This creates a back-action effect analogous to that 
producing a dynamical 'optical spring' in opto -mechanical 
systems 42 , and enables extremely steep trap barriers around z b 
(as in Fig. lc). In particular, the repulsive force is given by 



, cko e 

dz 



a ee (z) 



. where ^ 



is engineered to be large. However, 



the excited-state population a ee (z) ~ 4/^2 itself varies strongly 

with position through the spatially dependent detuning 

(for simplicity, here we omit the varying linewidth T(z)). 

Taking into account the spatial variation of cr ee around the 

equilibrium position, one finds an overall repulsive force that 
scales like (^) 2 , 



dCD e (z t )\ 2(z-Zt)(Tee(Zt) 



dz 



d a (z t ) 



(8) 



Results for trapping of caesium. Given the possible steepness of 
this barrier, an ideal limit for the overall potential is shown in 
Fig. 2a. The attractive part consists of the pure ground- state 
potential and arises from quantum fluctuations, while an infinite 
hard wall is created at z = z^ owing to the combination of strong 
resonant excited- state shifts (originating from the interaction of 
the atom with a classical image dipole) and opto -mechanical 
back-action. This idealized model enables one to understand the 



best possible scaling of the trap properties. In particular, the trap 
depth L/ depth (the energy required for a classical particle to 
become unbound) cannot exceed the value of the ground- state 
potential at z b , U m ^=h\Scj g (z b )\ (for example, (7 max ^ 60 mK 
when z b =10nm for a Caesium atom). This trapping potential 
can be easily quantized numerically. In Fig. 2b, we plot the 
quantum-binding energy U 0 (where 0 < U 0 < l/ depth < l/ max ) and 
position uncertainty Az for the motional ground- state wave- 
function. For the numerical results in Figs 2 and 3, we use atomic 
properties corresponding to Caesium (/L 0 = 852 nm, T 0 /2n = 5.2 
MHz, recoil frequency co r /r 0 = 4 x 10 ~ 4 ). The strength of the 
vacuum potentials is reflected in the strong confinement of the 
wavefunctions (for example, Az<4nm for distances z b <50nm), 
which are an order of magnitude smaller than in conventional 
optical traps. To contrast our results for Caesium with those of a 
less massive atom, we have considered Lithium as well. Results for 
Li similar to those for Cs in Fig. 2a,b show that at z b = 10 nm, the 
potential depth is smaller by roughly a factor of two, while the 
confinement Az sees roughly a fivefold increase. Generally, within 
the Wentzel-Kramers-Brillouin approximation and to lowest 
order in Az, the ground-state localization scales like k 0 Az~ 
((D r (k 0 z h ) 4 /r 0 ) 1/3 , while the binding energy U 0 &U max 
( 1 — 3 Az/z b ) can be nearly the entire classical potential depth. 

This ideal model describes qualitatively well the actual 
potentials over a large parameter regime. In Fig. 2c, we plot the 
ratio of the numerically obtained trap depth to the theoretical 
maximum, U dept h/U meLX , as functions of A p and dimensionless 
driving amplitude Q. Here, we have chosen a quality factor of 
Q=10 7 and a barrier position of z b =10nm. A trap depth 
comparable with (7 max is achievable over a large range of values. 
For low values of Q, the trap depth is weakened owing to escape 
over the potential barrier, while for large values, the trap becomes 
shallow as its position z t is pulled away from z b (as discussed 
previously). Two representative trap potentials, along with the 
ground- state wavefunctions and confinement Az are illustrated in 
Fig. 2d. 

Analysis of trap heating. We have identified four sources of 
motional heating— recoil heating, non-adiabatic motion of the 
atom within the trap, fluctuations of the atomic coherence and 
tunnelling over the finite-height potential barrier — that limit the 
trap lifetime in absence of external cooling mechanisms. Tun- 
nelling is exponentially suppressed with barrier height and thus is 
easily suppressed with increasing Rabi frequency. This calculation 
can be found in the Methods, while we discuss the more 
important mechanisms here. 

As in conventional optical trapping, the random momentum 
recoil imparted on the atom by scattered photons yields an 



. (M eff r 



R sc . Two 



increase in motional energy at the rate dE/dt 
qualitative differences emerge relative to free space, however. 
First, the photon- scattering rate R sc = T(r)<7 ee « i^z) 3 $ * s 

modified in the presence of dielectric surfaces. Second, the 
effective momentum M eff imparted by a photon is enhanced 
compared with the free-space momentum. This momentum 
scales at close distances like fc e ff ~ V^/z (see Methods), owing to 
emission of high-wavevector photons into the dielectric. Starting 
from the motional ground state, this heating causes the atom to 
become unbound over a timescale (k 0 z t ) 2 A p /(3(o 1 .). It should 
be noted that the heating rate dE/dt does not depend on the atom 
confinement (that is, whether we operate in the Lamb-Dicke 
limit), as the suppression of motional jumps in tight traps is 
compensated by the increase in energy gained per jump 43 . 
The increase in lifetime with trap position, r r oc z t 2 , is attributable 



entirely to the reduction in recoil energy 



(M eff ) 2 



associated with a 
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photon emitted into the surface. In particular, the decreasing trap 
depth L/depth oc z~ 3 versus distance is exactly compensated for by 
the reduced surface- enhanced emission rate r oc z t ~ 3 . 

Our derivation of the atomic potential assumes that the atomic 
coherence adjusts rapidly to the local trap properties, dcr ge /d£^0, 
which results in a purely conservative potential. As we now 
show, non-adiabatic corrections result in a momentum- 
dependent force, dp/dtccp (anti- damping). We can solve 
equation (6) for the atomic motion perturbatively by writing 

°ge(f) — °g?(z(0) + °geV)> where afj is the adiabatic solution. 
It depends only on the instantaneous position and is given by 
GgJ = ^jj^ in the weak saturation limit. Here, S c (z) = <5 a (z) + 
iT(z)/2 is the complex detuning. Substituting this solution 
into equation (6) yields a velocity- dependent correction to the 

coherence, <7gJ = (i3 c ) ~ 1 ^f- = — ^ This in turn yields a 
new contribution to the force, — ftp, where the anti-damping 
rate is given by ft « - « r <^^ [(r 2 - Ad\) f + 8<5 a rf ] 

evaluated at the trap position z t . Intuitively, the atom is anti- 
damped (/?>0) because the laser frequency is blue-detuned 
relative to the atomic frequency at z t , giving rise to preferential 
Stokes scattering over anti-Stokes. Aside from the term 
proportional to dI7dz, the heating rate is exactly analogous to 
that occurring in an opto -mechanical system in the blue-detuned 
regime 42 . We define a characteristic time T ad needed for the 
energy increase produced by anti-damping to equal the ground- 
state-binding energy U 0 . 

Thus far, we have focused on the motion of the atomic- dressed 
state, in which a trap emerges owing to weak mixing with the 
excited-state potential. Following a spontaneous emission event, 
however, the atom returns to the ground state and sees the pure 
ground-state potential for a transient time f t rans~ |^cfe)| ~ 1 until 
the atomic coherence is restored, during which the atomic wave 
packet accelerates toward the surface and gains energy. While a 
simulation of the full dynamics of the atom (including the spatial 
and internal degrees of freedom) is quite challenging, here we 
estimate the heating rate from these transient processes. In 
particular, we first find the spatial propagator under the pure 
ground-state potential, U(Zf, z b tf, t { ) (see Methods). Then, 
assuming that a spontaneous emission event occurs at t x — 0 with 
the atom initially in the motional ground state i/^fe), the new 
wavefunction that emerges once the coherence equilibrates is 
given by i//(Zf)~ fdzJJiZf, z b t tmnsy O^fe). The probability of 
atom loss P t rans per spontaneous emission event is estimated from 
the overlap between the initial and final states, P t rans — 1 — 
K^fki) | 2 > w ith a corresponding trap lifetime of T~ a ^ s = Ptrans^sc- 

The overall trap lifetime is obtained by adding the individual 
lifetimes in parallel, = (t~ 1 + t^ 1 + t"^ + t^), and is 
plotted in Fig. 2a. The overall trap lifetime exhibits a complicated 
dependence on Rabi frequency Q and detuning A p owing to the 
different scalings of the constituent heating mechanisms (see 
Fig. 3b). As an example, however, for external parameters 
A p = 2.5xl0 5 and Q/Q min = 1.5 x 10 4 for atomic Cs, one 
achieves a lifetime of T tota i^l5ms, ground-state confinement of 
Azi^lnm, trap depth of l/d ep th~ 10 mK and trap distance of 
z t ^15nm. This Rabi frequency corresponds to an incident 
intensity of just 0.05 mW um ~ 2 . To compare, nanofibre-based 
optical traps 8,11 employ guided intensities three orders of 
magnitude larger to achieve confinement and surface distances 
of Az~20nm (corresponding to trap frequencies of 
a? = 4nJ(Az) 2 ~ 100 kHz ) and £ t ~200nm, respectively. 
Trapping frequencies of co m /27itt40 MHz would be needed to 
reach single-nanometer confinement, far beyond what is realistic 
in optical traps. 



An additional heating mechanism, which does not appear in 
our simple model of an isotropic two-level atom, arises from the 
possibility of an excited state in a realistic multi-level atom to 
emit photons of different polarizations and transition into 
different states in a ground- state manifold. The different 
Clebsch-Gordan coefficients of these transitions, along with the 
symmetry breaking introduced by the nearby dielectric structure, 
imply that different excited states can experience different 
vacuum shifts. This would appear as different trapping potentials 
for the various dressed ground- state sub -levels absent of some 
special engineering. As a result, spontaneous changes in the sub- 
levels (such as arising from spontaneous photon scattering) 
would induce fluctuations of the trapping potential and lead to 
motional heating. One possibility to prevent these issues in an 
experiment is to utilize the effective two-level structure provided 
by a cycling transition. 

Discussion 

We have described a protocol for an atomic trap based on 
engineered vacuum forces, and have analysed in detail a model 
case of one- dimensional trapping near a Drude material, where 
analytical results are possible to obtain. The Drude response 
approximates well a number of metals, such as silver and gold, in 
the optical domain 44 . In practice, however, these metals are 
limited by their low-quality factors (Q<10 2 ), and the thermal 
fluctuations of such conducting materials can give rise to strong 
noise-induced trap heating and decoherence 45 . On the other hand, 
photonic crystal structures enable engineering of resonances 
through geometry, and quality factors approaching Q~10 7 have 
been observed 46 . We anticipate that our trapping protocol could 
be quite generally applied, and the resulting trap behaviour would 
qualitatively remain the same. In particular, the spatial 
dependence of the excited-state properties would no longer scale 
like (k 0 z) ~ 3 as in the case of a simple dielectric surface, but be 
characterized by a more complex spatial function /(r) that reflects 
the mode shape of the underlying structure 23 . However, the 
scaling of the excited-state shift (da> e ccQ/A p ) and emission 
(r oc Q/Ap) is generic to coupling between an atom and resonator 
of Lorentzian lineshape, and thus the capability of achieving 
highly stable traps with high-Q structures using our proposed 
techniques is maintained. Photonic crystals thus offer great 
flexibility in tailoring the properties of vacuum traps (such as 
dimensionality), which we plan to investigate in detail in future 
work. For example, it has been shown that dielectric gratings can 
produce quantum vacuum forces that trap atoms along the lateral 
directions of a dielectric interface 35 , which when combined with 
our previous techniques could yield full three-dimensional traps. 
We note that theoretical 23 and experimental efforts 13 have already 
been made to achieve realistic hybrid atom traps in photonic 
crystal structures, where a combination of vacuum forces and 
optical forces is required to stabilize the trap in all dimensions. In 
addition, while we have focused on trapping of atoms weakly 
dressed by an external laser, it would be interesting to investigate 
the feasibility of laser-free traps, where atoms pumped into 
metastable electronic levels are stably confined using engineered 
resonant shifts alone. 

The ability to create atomic traps with parameter sets (such as 
depth, confinement and proximity to surfaces) that are not 
possible in conventional traps (whether optical, electrical or 
magnetic) should lead to a number of intriguing applications. For 
example, it has been proposed that the trapping of atoms near 
dielectric surfaces with nanoscale features (such as a sub- 
wavelength lattice) would facilitate large interaction strengths 
and long-range interactions for ultracold atoms 39 . However, the 
major limitation of optical trapping is the divergent intensity 
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needed to overcome the ground- state vacuum forces of these 
structures. The use of resonantly enhanced excited-state repulsion 
to create a trap could significantly reduce the power requirements 
and enable smaller lattice constants. Furthermore, traps with 
localization on the nanometer scale can induce inter- atomic forces 
that are of comparable strength with molecular van der Waals 
forces, which enables novel opportunities to control ultracold 
atomic collisions 20 and realize exotic interactions such as p-wave 
scattering of fermions 47 . Atoms trapped near surfaces could also 
act as exquisite nanoscale probes of surface physics, including the 
precision measurement of vacuum forces themselves. 

Methods 

Electromagnetic field quantization. In this section, we briefly describe the 
quantization scheme of the electromagnetic field in the presence of arbitrary linear 
dielectric media, which follows that of ref. 27. The free field is characterized 
by a set of bosonic modes with frequencies co and annihilation operators aj (r, co) 
(j = x, y, z), with a corresponding Hamiltonian 

Hf = hY^ J dr J dwcoaj (r,(o)aj(r,a>). (9) 

, The bosonic operators satisfy canonical commutation relations [a ; -(r,co), 
a[ (Y, co')] = djkd{r — r')d(co — co') and are physically associated with noise 
polarization sources in the material. The electric field operator at frequency co is 
driven by the sources at the same frequency and is given in the Coulomb gauge by 27 

Ej(r, co) = ip 0 co 2 J dr'Gjk(r, r', co) ^-^-lme(r r , coja^r', co). (10) 

Here and in the following, it is assumed that all repeated vector or tensor 
indices are summed over, for example, GyOj = J2^ij^j- c(r,co) is the frequency- 

; 

dependent permittivity at position r, and Gy is the dyadic Green's function which is 
the (gauge-invariant) solution to 



(VxVx)- ^£(r,a>) 



G(r,r',<L>) = 5(r-r') <g> I. 



(11) 



The total electric field is E(r) = fdco E(r, co) + h.c. Similarly, the magnetic field 
is B(r) = f dco ± Vx£(r, co) + h.c. It can be verified that these expressions 
preserve the same field commutation relations as in free- space field quantization. 

We now consider a two-level atom at position r a and with ground and excited 
states \g}, \e}, respectively, interacting with the electromagnetic field within the 
electric dipole approximation. The interaction Hamiltonian is given by 



-d-E(r a ) 



where p { is the atomic dipole matrix element. 



(12) 
(13) 



Plane-wave expansion. We consider an interface with vacuum in the region z>0 
and a dielectric with permittivity e{co) in the region z<0. The solution to 
equation (11) in the region z>0 can be written as the sum of a free component 
(that is, the solution in isotropic vacuum) and a scattered component. In the region 
0<z<z' y they can be expanded in plane waves, which yields 



Gfre 



— f 

An 2 co 2 J 



5 (z - z!)zz + ^- e ik± {z ' z) k 0 x k 0 x 



T 2 C0 2 j ( 



Jc 2 

= 87T 2 , 



r p (k\\z-k±k\\ )(k\\z + k±k\\ )+r s — (zxk\\ )(zxk\\) 



(14) 



(15) 



Here k|| = (k x , k y ) and p = (x, y) are the wavevector and position along the 
direction parallel to the interface, while the full wavevector is given by k 0 = 

k|| — kj_z with the relation k\ J rk\ = (co/c) 2 . Further defining q± = 



e(co)(co/c) — k 2 \ as the perpendicular component of the wavevector in the 

sdium, r s = \ L ~ q± and r v = \ L<Z ~ L 
and p polarized waves, respectively. 



l ±e - q± are the Fresnel reflection coefficients for 5 



Ground-state interactions. The interaction Hamiltonian H a{ couples the ground 
state of the bare system \g, 0 > (consisting of the atomic ground state and vacuum) 

to states \ e, l T ,<aj) = ( r ? w ) l e >0)- Within second-order perturbation theory, 



this interaction induces a shift of the bare ground- state energy by an amount 



5o) g (r a ) = - ^ J2 J dr / 



dco 



(e,l r , m ,;|H af |g,0)| 2 



CO + COq 



(16) 



where co 0 is the bare resonance frequency of the atom. Substituting equations (10) 
and (13) into the expression above, and using analyticity of the integrand to rotate 
the frequency integral onto the positive imaginary axis co = iu, one finds 27 

oo 

/2 
du r^T 2 & i G sc,y ( r a , r a , iu) pj . (17) 
( J ) q ~~ | hi 

0 

Here, we have only included the scattered component of the Green's function, 
as the free component gives a contribution that is independent of position. While 
we have considered a two-level atom for simplicity, a more realistic model of an 
isotropic atom consists of multiple excited states and equal polarizabilities in each 
direction (that is, the ground state can emit virtual photons of any polarization 
with equal strength). Since the contribution to the ground- state shift from each 
excited-state transition is additive, equation (17) can be modified to account for 
this by setting p { = p 0 for all i and associating colpl/(37ie Q hc 3 ) = r 0 with the 
free- space spontaneous emission rate of any excited state. This yields equation (1). 

Excited-state interactions. We want to analyse the effective excited-state prop- 
erties in the specific case that an atom is weakly driven by a classical field of 
frequency co L . We thus transform to a rotating frame where the free atomic evo- 
lution is given by H a = — 3a ee and 3 = co L — co 0 is the detuning between the laser 
frequency and atomic resonance frequency co 0 . 

The standard technique to derive the excited- state shift and decay rate follows 
that of the ground state and uses time-independent second-order perturbation 
theory 27 . In this approach, however, one cannot find the jump operator associated 
with photon emission, which is needed to calculate the effect of emission on atomic 
motion (for example, recoil heating). To rectify this, we employ an alternative 
approach based on deriving an atomic master equation within the Born-Markov 
approximation 48 . In particular, the equation of motion for the reduced atomic 
density matrix p a is given by 



Pa : 



x 

f /d 4 



Tr t / dT[H af (f), [Haf(<-T),p a ® 1 0> <0 1 ] 



(18) 



where H a f is the interaction Hamiltonian in the interaction picture (with respect to 
the free Hamiltonian H 0 = H a + H { ) and |0> is the electromagnetic vacuum state. 
Any exponential of time appearing in the integral can be evaluated using the 
relation dx e wn = nd(co) + iV\ Q . Using the field expansion of equation (10) in 
H a f and following some manipulation, equation (18) can be written in the form 



-[fW a ] + £[p a ]. 



(19) 



The effective Hamiltonian describes coherent interactions between the atom 
and vacuum field and takes the form 



H e ff = fi(5(L> e (r a )<7 e£ 



(20) 



which can be interpreted as a position-dependent energy shift of the excited state. 
Taking a simple isotropic two -level model for an atom, the excited- state shift is 
found to be 



Sco e (r a ) - 



cT 0 



x 

/ 



du 



03>l+U 2 



Tr G sc (r a , r a , iu) - 



TqUC 



Tr Re G sc (r a , r a , co L ), 
(21) 



as stated in the main text. 

The Liouvillian term C[p a ] can be written in a Lindblad form 



£[p a ]= W dk0 k;Pa 4-^(4^+Pa4^)- 

; J 

The jump operators take the form 



(22) 



fi J ~(2nf 2 

The excited- state spontaneous emission rate is in turn given through the 
relation / 0* O kj = hr(r a )cr ee , or 



r(r.) = 



-Tr ImG sc (r a ,r a ,co L ). 



(23) 



(24) 



This again averages over dipole orientations and reproduces equation (3), with 
the replacement of co 0 -> co L . To derive this, we have used the identity 



oj 2 f 

\ dr // Ime(r // ,a;)Gy(r,r // ,a;)G;(r // ,r / ,ft;) = 



: Im G/fc(r, r', co)- 



(25) 



The expressions for the excited- state emission rate T(r a ) and level shift dco e (r a ) 
agree with previous derivations 27 , if the laser frequency co L is replaced by the 
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atomic resonance frequency co 0 . The appearance of co L physically describes the 
effect of the weakly driven atom (that is, a linear optical element) interacting with 
its own Rayleigh- scattered field, which can be altered owing to the presence of 
nearby dielectric surfaces. This is in accordance with the well-known result for the 
fluorescence spectrum of a driven atom, which is dominated by the Rayleigh 
contribution for atomic populations a ee « 1 ( re ^- 4 8). More importantly, we have 
obtained for the first time a decomposition of the jump operator associated with 
photon emission for arbitrary dielectric surroundings, which enables one to 
calculate recoil heating. 

Short-distance limit. We now derive the short- distance expansion of the excited- 
state properties. The resonant term in equations (2) and (3) depend on the quantity 
Tr G sc (r a , r a , cl> l ), which for planar interfaces is given by 



oo 

TrG sc (r a ,r a> co L ) =^ J dk H ^-e 2ik ^ (r, + ^ (k\ 



(26) 



For short distances, the expression is dominated by large k\\&ik±, which gives 



rise to a long exponential tail e 2lk±z 
k\\, this tail contributes as 



. Expanding to second order for large 



TrG sc (r a ,r a ,a> L ) « — 
In 



X 



dk\\ e 



c 2 k\ e -l 
co 2 e + 1 



e(e-l) 



(27) 



We now apply these results to the Drude model. First, we note the exact relation 



e(co): 



2Q 2 

(Q + A p 7+7(Q + A p ) 



. For frequencies within a range A p < Q of the plasmon 



resonance, one readily finds that the response |— j « — f ^ A is well approximated 
by a complex Lorentzian. For larger detunings A p > Q, the exact expression of e(co) 
must be used. Furthermore, within the range A p < Q, one must compare the size of 
the first and second terms on the right-hand side of equation (27). The first term 
dominates the integral provided that k 0 z<^A p /Q, in which case 

TrG sc (r a ,r a ,o; L ) « - ^ 3 ^ , (28) 

which reproduces the asymptotic expressions for the excited- state shift and 
linewidth in the main text. 

Tunnelling. The atomic escape rate from the metastable potential formed by the 
vacuum trapping scheme owing to quantum tunnelling is calculated in the 
Wentzel-Kramers-Brillouin approximation in two steps. First, we find an 
approximation to the real part of the bound- state energy tico h . This is implemented 
by replacing the actual metastable-dressed-state potential U d (z) with a stable one, 
U s (z), according to the formula (also see Fig. 4) 

17. (z) = U d (z) (z>z max ), (29) 



N 



^d(^max) — ^max)' ( 


^rna 


X 




L 



Figure 4 | Tunnelling. Black curve: metastable-dressed-state potential 
(J d (z) (in arbitrary units) seen by the atom. The maximum of the 
metastable barrier is located at z max . Red curve: stable potential (J s (z) used 
to approximately calculate the binding energy frco b (blue line) of the 
atom in the trap. The classical turning points z lr 2 of a particle with this 
energy, which are used to evaluate the tunnelling rate in the Wentzel- 
Kramers-Brillouin approximation, are labelled as well. 



Here, z max is the position corresponding to the maximum of the metastable 
barrier. The ground-state energy of U s (z) is solved numerically. 

Once the approximate ground- state-binding energy is obtained, the tunnelling 
probability P t for a particle hitting the metastable barrier is calculated using the 
Wentzel-Kramers-Brillouin approximation, 



P t = 4 exp 



2m(Ud(z) -hoj h/ 



(31) 



Here, z lj2 are the solutions to U d (z) = hco h and are thus the classical turning 
points of a particle with energy hco h . The rate that the trapped atom collides with 
the barrier is approximated by ApJmAz, where the position and momentum 
uncertainties are obtained from the numerical ground- state solution. We then 
estimate the inverse of the tunnelling-limited lifetime to be t~ 



L tunnel mAz ' 



Recoil heating. Recoil heating can be calculated by treating the atomic position r a 
in equation (23) as an operator. For notational simplicity, we shall consider recoil 
heating only along the trapping direction z, although our results can be easily 
generalized. The increase in momentum uncertainty owing to photon scattering is 
given by 



dt 



(Pz>=Tr(p^ a ). 



(32) 



Here we will focus solely on the Liouvillian term in the density matrix evolution 
(that is, momentum incurred from spontaneous emission). It can be shown that the 
Liouvillian term produces no net force, that is, ^ (p z ) =0. Using the fact that 
[pz,/( z a)] = — Jjr for any function/, one can derive the following general 
expression for the recoil heating rate in terms of Green's functions, 



dt 



(pi) = H^cdIp^j (d\ + 2d,d 2 - d 2 2 ) (Im Gy (z a , z ai cn L )a ee ). (33) 

Here, the derivatives d\ and d 2 act on the first and second spatial arguments of 
the Green's functions, respectively, and we have suppressed the atomic spatial 
variables in the directions that are not of interest. We further assume that the 
internal and spatial degrees of freedom can be de-correlated, (Gy<7 ee ) ~ (Gy)a ee . 
This is valid in our regime of interest where properties such as the excited- state 
decay rate have negligible variation over the spatial extent of the atomic wave 
packets. As in the main text, we use the notation ay = (ay) to denote the 
expectation value of atomic operators. Formally, we can write the above equation in 
the convenient form 



dt 



(Pi 



(fifceff) r(r a )(7 e£ 



(34) 



This form is intuitive as it describes the momentum increase arising from a 
random walk process, where the atom experiences random momentum kicks of 
size ± M e ff at a rate F(r a )a ee (that is, the photon- scattering rate). fc e ff is a derived 
quantity that thus characterizes the effective momentum associated with scattered 
photons in the vicinity of a dielectric surface. 

As a simple example, we consider an atom in free space polarized along x. 
Evaluating equation (33) wit h the free-space Green's function Gf ree; xx from 
equation (14) yields fc e ff = y/2/5(co^/c) (and T(r a ) = r 0 ), recovering the known 
result 43 . This recoil momentum is smaller than the full momentum co L /c of the 
scattered photon because the atom emits in a dipole pattern, but only the 
projection of the photon momentum along z contributes to motional heating in 
that direction. Applying the same formalism to an atom near a dielectric described 
by the Drude model, we find that the scattered component of the Green's function 
dominates the recoil heating, leading to the asymptotic expansion of k t & w 
given in the main text. In all of our numerical calculations, the full Green's function 
is used rather than asymptotic results. 

Transient heating. Following a spontaneous emission event, the spatial wave- 
function of the atom is temporarily untrapped and evolves under the pure ground- 
state potential for a characteristic time £ trans before the internal dynamics equili- 
brates. Here we derive an approximate propagator U(zf, z b f f , for the wave- 
function evolving under the ground- state potential. 

The problem significantly simplifies if the ground- state potential U g (z) = ftco s (z) 
is linearized around the atom trapping position z = z t , 



^ n +U s (z t ) + U' g (z t )(z-z t ). 



(35) 



This linearization is justified in our regime of interest where the atom is tightly 
trapped and the wavefunction undergoes little evolution in the time t trans . The 
constant energy and position offsets in the Hamiltonian play no role, and thus we 
can consider the simplified Hamiltonian H = ^ + U' g (z t )z. The propagator for 
this Hamiltonian is given by 



U(Zf, Z{, t, 0) = Ufree(Zf,Zi, t) exp 



iU't(zf + Zi) i(urft : 



in 



24hm 



(36) 



where Uf ae (zf, 
free particle 49 . 



is the well-known propagator for a 
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In our simulations, the ground- state wavefunction for an arbitrary trapping 
potential is obtained numerically on a lattice, and so the propagator must be 
spatially discretized as well. In particular, while the exponential term in 
equation (36) is well-behaved on a lattice, the free propagation l/ free is not at short 
times. We instead replace l/ free with the propagator D mn (t) for the free discrete 
Schrodinger equation, 



dil/ n 



■■m n + l(t)-^n(t)+*n-M), 



(37) 



where /? = hllma 2 , n is the lattice site index, and a is the lattice constant. The 
propagator for this equation (denned as the matrix satisfying i// m (f) = D mn \J/ n (0)) is 
given by 



D mn (t) = i m - n U-n(2pt)e 



- apt 



(38) 



where / n is the ft-th order Bessel function. 
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